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ABSTRACT 

Aims. We calculate parameters A and B of the Baym-Pines model of the hydro-elastic equilibrium of rotating neutron 
stars. Parameter A determines the energy increase of a non-rotating star due to a quadrupolar deformation of its shape. 
Parameter B determines residual quadrupolar deformation due to the crustal shear strain in a neutron star that spun 
down to a non-rotating state. 

Methods. The calculations of A axe based on precise numerical 2-D calculations for rotating neutron stars with the 
reahstic equations of state (EOSs) of dense matter. An approximate, but quite precise, formula for B is used, which 
allows us to separate the contribution of the crust from the dependence on the stellar mass AI and radius R. The elastic 
shear strain distribution within the crust is modeled following Cutler et al. (2003). Realistic EOSs of neutron star cores 
are used, some of them with high-density softening due to the appearance of hyperons or a phase transition to an exotic 
state. 

Results. The values A{M) and B{M) were calculated for 0.2 Mq < M < 0.9 Mmax (where Mmax is the maximum 
allowable mass) for seven EOSs of neutron star core, combined with several crust models. A standard formula based on 
the incompressible fluid model is shown to severely underestimate the value of A. For M < 0.7 Mq the values of A{M) 
are nearly EOS-independent and are given (within a few percent) by a universal formula A — 3.87 {M / MqY x 10^^ erg. 
We derive the scaling of B with respect to R and M, also valid for a thick crust. We show that B for accreted crust 
strongly depends on pycnonuclear fusions at p > 10^^ g cm~'^. 

Key words, dense matter - equation of state - stars: neutron 



' 1. Introduction 

The solid crust of neutron stars forms quickly after their 
, birth. It plays an important role in neutron-star dynam- 
' ics. Apart from pressure, the stresses in the crust are pro- 
\ duced by an elastic shear strain; therefore, the crust can 
■ have a non-axial deformation of mass distribution, which al- 
' lows a rotating neutron star to radiate gravitational waves. 
, Moreover, the inner neutron-star crust, with a crystal lat- 
' tice of nuclei immersed in neutron superfluid, can be a 
site of pinning the neutron vortices to nuclei. A massive 
breaking of this pinning after a maximal strain has been 
reached is thought to be responsible for triggering the pul- 
sar glitches. Finally, the presence of the crust makes a free 
precession of a rotating neutron star possible. 

The hydrostatics of the neutron star crust is more 
complicated than for an ideal fluid, because of the elas- 
tic strain that produces an anisotropic contribution to 
the stress tensor. To find the mechanical equilibria of a 
crust with non- vanishing shear strain, one has to solve par- 
tial differential equations of the hydro-elastic equilibrium, 
with appropriate boundary conditions on the stellar sur- 
face and on the crust - liquid core boundary. A simple 
"one-parameter model" of a neutron-star structure with 
crustal elast ic strain and it s appl ication to pulsar glitches 
was given bv lBavm fc Pinej (|l97lh . The Baym-Pines model 
was reformulated in the framework of General Relativity 
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by ICarter fc QuintanBl (|l975l). A neutron -star precession 
model was developed bv iPines fc ShahamI (1972a, b), who 
used methods originally developed to d escribe rotation 
of the Earth (|Munk fc McDonaldl Il960l and references 
therein). In these models, the neutron star is treated as 
a solid shell (crust) enveloping a liquid core. More recently, 
updated revi ews of models of processin g neutron stars was 
presented bv I Jones fc AnderssonI (|200lD . The crustal rigid- 
ity, resulting from elastic strain in the deformed crust 
of neutron stars, was r ecently reexamined, using realistic 
dense matter models bv lCutler et al.l (|2003[ ). 

The change in the stellar energy due to star's quadrupo- 
lar deformation is d etermined by the parameter A of 
iBavm fc Pines! (|1971h . In the present paper we calculate 
A versus stellar mass M, using quadratic expansions of 
exact quantities obtained in the precise 2-D simulations. 
To this aim, we have developed a systematic procedure 
enabling us to calculate a quadratic response of neutron 
star structure to a rigid rotation. Using a set of EOSs, we 
study the EOSs dependence oiA{M). Then we calculate the 
parameter B measuring the elastic strain contribution to 
the hydro-elastic equilibrium of a rotating star with crust. 
In doing this, we avoid solving the complete hydro-elastic 
equilibrium e quati ons, by relying on results obtained by 
ICutler et al.l ()2003[ ). and adapting them to different EOSs 
of the crust and core. We then study the scaling of B with 
the mass and the radius of neutron star and investigate the 
differences between the ground-state and accreted crusts. 
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The plan of the paper is as follows. The model of Baym 
and Pines is briefly described in Sect. [21 where the param- 
eters A and B are defined. Calculation of A, based on the 
2-D simulations for rigidly rotating neutron stars with re- 
alistic EOSs, is described in Sect. [3] The shear modulus 
of the crusts is briefly reviewed in Sect. |4l where we also 
discuss the differences between the ground-state and ac- 
creted crusts. An approximation suitable for calculating the 
crustal rigidity parameter, valid also beyond the limit of a 
thin crust, is proposed in Sect. El Existing calculations of 
B, and our results for this quantity, are described in Sects. 
[6l[7l Scaling of B with M and R is briefly discussed in Sect. 
[H Numerical results for crustal rigidity versus stellar mass 
are presented in Sect. [HI Discussion of our results and con- 
clusions are presented in Sect. [TOl 



2. The Baym-Pines model 

We restrict ourselves to neutron star stars with an ellip- 
soidal (quadrupole) deformation of an axially symmetric 
crust. The rest mass of the star is fixed and equal to Mb. 



the dominant, ideal liquid term cr. 



liq 



The stress tensor of dense matter is — cr'''^ -I- crfi^, where 

J IJ LI' 



-P5ij, P is the pres- 



sure, and an elastic shear stress component afj' (see Sect. 
5) exists only within the solid crust. Correspondingly, the 
energy of a configuration \s E = i?'"^ -I- E^"^. 

Non-sphericity of an axially symmetric rotating con- 
figuration results in a moment of inertia tensor 7^ = 
diag{Ii, Ii, I3). The quadrupolar deformation is fully de- 
scribed by the oblateness parameter 



/3-/0 

and diagonal elements of lij are 

/i =/2 =/o(l-e/2) , /3=/o(l + e), 



(1) 



(2) 



so that the trace of lij (invariant with respect to rotations 
of the coordinate system) is S/q. 

Our definition of oblateness parameter requires a com- 
ment. It is the same as us ed bv lCutler et all (I2003D : how- 
ever, it differs from that of lBavrn &: Pinej l l971 ). who de- 
fined it in terms of the crust moment of inertia, cbp = 
(/c3 — Ico)/IcO- For stars made of incompressible fluid, 
e = eBp; however, for realistic stars, made of compressible 
matter, e > e^p- 

Two reference configurations are defined as follows. 

Co - non-rotating and relaxed. It is spherical, with 
E^^{Co) = (no shear stress, completely relaxed) and has 
energy Eq. Its moment of inertia tensor is X*^- = laSij. 

C - rotating and relaxed. It is rigidly rotating, at rotation 
frequency Cl = 27r/ (rotation period). This configuration is 
also completely relaxed, so that E^^{C) = 0. 

An example of a reference configuration C is a newly 
born neutron star that is already rigidly rotating but still 
completely liquid. Then the star spins down, and its elliptic- 
ity decreases to e < e. However, the crust solidifies quickly, 
and elastic strain accumulates within the crust, which is 
opposing the decrease in e. 



Let us consider a rotating configuration C, with rotation 
frequency fl, and born as C. Assume that there is a continu- 
ous spin-down transformation C — > C (no crust breaking). 
Then, a "residua l oblateness" , which adds to that resulting 
from rotation, is (Bavm fc Pinesl . ri97ll : ICarter fc Quintanal 
,1975, : .Cutler et al.. . ^03) 

e, = be . (3) 

Here, 6 is a "rigidit y parameter" (called so by 
iJones fc AnderssonI 1200 1[ ) , which depends on the elastic 
properties of the crust, on the EOS of neutron star mat- 
ter, and on the mass M of the reference configuration Cq. 
The name "residual oblateness" stems from this being the 
ellipticity that remains after a complete spin down toil = 0. 

Within the Baym-Pines model, the value of b can be 
calculated by considering the limit of a non-rotating star. 
The quadratic approximation to the stellar energy is then 



E — Eq 



Ae^ 



B{e 



(4) 



where parameter A is calculated from the quadratic 
response of a non-rotating neutron star energy to a 
quadrupole deformation. The last term on the righthand- 
side is E^^, verifying E^^iC) = 0. The equilibrium of a 
non-rotating star is reached at minimum E, which implies 

B B 

In their calculation of b, ICutler et al.l (|2003[ ) rely directly 
on the difference Al = I3 — Iq = Xae, calculated in two 
cases. The value AIq = Joe is obtained for a fully relaxed 
configuration rotating at fi, while the residual A/^ = /oCr 
results from a continuous spin-down from il to zero. Within 
the quadratic approximation in fl, one gets A/d/A/o = b = 

Early calculations for a 1.4 Mq n eutron s tar hav e given 
A ~ 10^^ erg and B ^ 10*^ erg ()Pandharipande et al.l 
1X9761 ). so that b ~ B/A ~ 10~^. Recent calculations for a 
1.4 Mq neutron star with realistic EOSs of neutron star 
crust and core yields a in uch lower value 6 = (1.7 — 2.5) x 
10-'^ (jCutler et all . 120031 ). D ependence of A and B on the 
stellar mass was studied by ICarlini fc TrevesI (Il989l). As 
we wi ll see, the approximations used by ICarlini fc TrevesI 
were rather unrealistic. 

3. Rotational deformation of a relaxed 
configuration 

We intend to calculate parameter ^ as a function of the 
gravitational mass of the non-rotating configuration of the 
same baryon mass, Mb. Consider configuration rotating 
rigidly at angular frequency f2, as measured by a distant ob- 
server. In what follows, we use the results of exact 2-D cal- 
culations. Einstein equations for stationary axi-symmetric 
space-time are solved using the rotstar code, a part of 
the LDRENE public-domain C-f + scientific library based on 
spectral methods of solution of partial differential equations 
(see http://www.lorene.obspm.fr/). The stellar angular 
momentum is J and mass-energy is Mc^. In what follows, 
in the spirit of the Baym-Pines model, we restrict ourselves 
to a slow rotation limit; i.e., we keep only lowest powers of 
fl. To extract relevant coefficients, we have chosen a non- 
dimensional expansion parameter equal to the angular 
frequency fl divided by the angular frequency for a star 
rotating at 10ms period, = rj/ilioms = ri/(27r • 100 s~^). 
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Table 1. Equations of state of neutron star core. 



EOS 


model 


reference 


BPAL12 


N energy density functional 


a 


FPS 


N energy density functional 


b 


GN3 


N relativistic mean field 


c 


DH 


N energy density functional 


d 


APR 


N variational theory 


e 


BGNIHI 


NH, energy density functional 


f 


GNHQm2 


NH + mixed baryon-quark state 


g 



NH - hyperons, nucleons, and 
[Bornbacil 



Note: N - nucleons and leptons. 
leptons. 

References for the EOS : a - iBombacil (Il995l) ; b 
iPandharipande fc Ravenha.111 lll989f ): c - iGlendennind (119851 ): 
d - iDouchin fc Haensell (l2001h: e - A18J +UI X* model of 
lAkmal et all l| 19981 ): f - iBalberg fc Gall ||1997D : g - iGlendennirig 
l|2000l ). 



3.1. Coefficients a of the star's response to rotation 





0.5 1 1.5 3 

M [MJ 



Fig. 1. Parameters as, ctj, and ad versus stellar mass M, 
calculated for the DH EOS (Table 1). For further explana- 
tions see the text. 

The non-sphericity of the mass-distribution is measured 
by ellipticity e. We calculate and expand I3 using 

h ^ ^ = Io{l + ain^) , (6) 

which yields 

e = aiTf . (7) 

Simultaneously, we can parametrize the geometrical oblate- 
ness of the shape of rotating star by 

2 

I - rpoi/r^q^ asfl , (8) 

where rpoi (j'cq) is the polar (equatorial) radial coordinate 
of rotating star. Let us now consider the decomposition 



Fig. 2. Same as Fig. [T] but calculated for the BGNIHI 
EOS (Table 1). 

of the total mass-energy Mc^ into the static part and ki- 
netic energy. The kine tic energy is defined b y analogy to 
the Newtonian theory (|Friedman et al.l Il986l) 

^^1^^= l^^' ■ (9) 

Therefore, 

Mc^ = M^°'>c^ + Edcf + T , (10) 

where M^o^c^ is mass-energy of non-rotating configuration, 
and i?dcf is the increase in intrinsic energy of the star (i.e., 
calculated in the star's reference system), as compared to 
the non-rotating (static) case. To the lowest order, the lat- 
ter quantity can be parametrized as 

-Bdct = ad (1 - r-poi/rcq)^ . (11) 

The parameters a were calculated, using exact 2-D mod- 
els of rotating neutron stars, for several realistic EOSs 
of neutron star cores. These EOSs are listed in Table 1. 
As the EOS of neu t ron st ar crust, we used the model of 
iDouchin fc Haen^el (|200lh (hereafter referred to as DH). 
Examples of the dependence of the a parameters on the 
gravitational mass of neutron star are shown in Figs. [1] (DH 
EOS) and El (BGNIHI EOS). The DH EOS assumes that 
the matter is composed exclusively of nucleons and leptons, 
and it leads to a monotonous behavior of all a-s. Namely, 
a I and decrease with increasing M. Such behavior is 
quite natural, because as and a/ measure the response of 
the shape and mass distribution to a rotation at a fixed fl. 
This response (polar fiattening, anisotropy of the moment 
of inertia tensor) weakens with increasing mass, because 
gravitational pull, which tends to concentrate the matter 
towards the center, becomes stronger. On the contrary, the 
increase in energy, relative to the spherical shape, due to a 
given ellipticity of the star, is measured by ad- This energy 
rapidly increases with star mass: gravitational binding of 



4 



J. L. Zdunik et al.: Deformation and crustal rigidity of rotating neutron stars 



the star rises with M, and the star becomes stiffer to polar 
flattening. 

There are some bounds on M resulting from the 
high-precision constraint imposed on the 2-D simulations. 
Calculations of as had to be stopped at M ~ 0.9 Mmax, 
because for M « Mmax, the closeness of relativistic insta- 

— 2 

bility makes precise 2-D calculation of linear response to VL 
very difficult. On the other hand, M could not be too low, 
because for low M the rotating star becomes very "soft" 
to rotational deformation, which is reflected in very high 
values of as and a^ . All in all, very precise calculation of 
a-s could be done only for 0.2 Mq < M < 0.9 Afmax-Q 

The BGNIHI EOS takes the high-density softening due 
to the appearance of hyperons into account. However, for 
M ^ 1.35 Mq stellar core does not contain hyperons, 
and the a{M) dependence is monotonous, similar to the 
DH EOS. For M > 1.35 Af© hyperons soften the neutron 
star core, resulting in a specific M{R) dependence. Namely, 
M increases very weakly with a decreasing R (flat M{R) 
curve). In contrast, one notes a strong increase in the bind- 
ing energy of the star, and a decrease in the moment of 
inertia with increasing M, Fig. [31 The effects of the hy- 
peron softening on a{M)-s are therefore the following. The 
function as (aa) decreases (increases) more steeply with M, 
because the star becomes more compact and its binding in- 
creases rapidly with M . The most characteristic, however, 
is the effect on ai{M), which increases rapidly by a factor 
of two, reaching a maximum at M ~ 1.5 Mq, and then falls 
monotonously. This specific behavior is to some extent due 
to a rapid drop in / for ill ^ 1.4 Mq. 



3.2. Energy and oblateness - coefficient A 

Using Eqs. ([7]), ([8]), and (fTTjl we obtain final expression for 
A in terms of three as: 



A^ ad [as/aiY 



(12) 




2.5 



2 - 



E 1-5 
o 

BO 

o 

^ 1 



0.5 



~I [ [ T" 



APR / 




GMHl 



M [MJ 



Fig. 3. Moment of inertia for slow, rigid rotation, I — Iq, 
versus stehar mass, M, for the APR and BGNIHI EOSl 
(Table 1). 



Fig. 4. Parameter A versus stellar mass M for the APR, 
DH, and GN3 EOSs (Table 1). For more explanations see 
the text. 
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FPS / 




BPAL12 / 




/V''\ BGNIHI - 




/ A/ I 




/// \ < 




/ // , / 




/ // .. ' ; 




/'/ ■■ ' 








/y ■■ 




GNnQm2 V 
1 < < < < 1 < < < < 1 < , , 



0.5 



1.5 



M [MJ 



^ Let us remind that non- rotating cold neutron stars Mn 
0.1 Mq (|HaenseI et al.ll2002D 



Fig. 5. Same as Fig.Hbut for the BPAL12, FPS, BGNIHI, 
and GNHQm2 EOSs (Table 1). For more explanations see 
the text. 
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Plots of A{M) for seven realistic EOSs from Table 1 are 
shown in Figs. [3] and O As already explained in Sect. 13. ll 
high precision 2-D calculations were only possible within a 
limited range of stellar mass, 0.2 Mq ^ M ^ 0.9 M„iax- 
Notice that EOSs in Fig. 3] have no high-density soften- 
ing. They lead to monotonously increasing A{M), with a 
maximum scatter at given M characterized by a factor of 
two. 

In Fig. [5] we compared A{M) for two EOSs with high- 
density softening, due to hyperonization (BGNIHI) and to 
a transition to a mixed hadron-quark phase (GNHQm2), 
with two EOSs without such a softening. Hyperon soft- 
ening implies a dramatic drop (by a factor of five!) for 
M ^ 1.4 Mq. The appearance of a mixed-phase core 
leads to quasi-constancy of A for BNHQm2 EOS near 

M ~ 1 Mq. 

Up to now, A{M) has been calculated using the model of 
incompr essible-fluid stars in t h e Newtonian theory of grav- 
itation dBavm fc Pined llQTlt iPandharipande et all 119761: 



ICarlini fc TrevesI 119891 : 1 Jo nes fc Anderssonll2001h . For such 
a model, th e formula for A{M) was originally derived by 
ILo^ (Il920l) . 

As we see in Fig. [6l this formula badly underestimates 
A{M) for realistic neutron stars with a high-density soft- 
ening, and it is unable to reproduce a strong drop, due to 
a high-density softening of the EOS. 



25 



20 - 



15 



10 



APR 




BGN1 



M [MJ 

Fig. 6. The coefficient A versus stellar mass, M, for the 
APR (thick dash line) and BGNIHI (thick sohd hue) EOSs 
(Table 1). Thin hues (dash - APR, solid - BGNIHI) cor- 
respond to a standard incompressible-fluid approximation, 
Eq. (ini). 



For a sufficiently low stellar mass, one can describe the 
A{M) relations by a universal formula, because the influ- 
ence of the dense-matter EOS becomes less and less impor- 
tant with decreasing M. Consequently, the A{M) relations 



for different EOSs converge to a single one. An approximate 
formula (fltted in the 0.1 — 0.7 Mq range) reads 



A{M) = ai {M/MqP 10^'^ erg 



(14) 



with ai = 3.87 and bi = 7/3. The fltting formula works 
very well, with accuracy of the order of one per cent, for 
masses less than 0.7Mo. Its quality worsens for higher 
masses when the A{M) relations become EOS-dependent. 
Similarly, one can approximate the numbers given by the 
Newtonian Love's expression, where the true values of M 
and R are used, by a universal function Alovc(-M) analo- 
gous to Eq. (fH)) . but with ai = 2.68 and bi — 2. Such a 
formula for Ai^ovc{M) is quite precise (within a few percent) 
for masses lower than O.7M0. 

4. Shear modulus of the crust 



We follow the notations of iHaensel et al.l (j2007[ ). 
Deformation of the crust with respect to the equilib- 
rium (ground state) configuration of ions (nuclei) in the 
lattice implies displacement of ions into new positions 
r' = r + u, where u — u{r) is the displacement vector 
(field). In the macroscopic (continuum) limit, r and u 
are treated as continuous fields. Nonvanishing u induces 
deformation energy density S^d — — So, where Sq 

corresponds to the non-deformed (equilibrium) crust. 
Non-zero u induces forces that tend to return a matter 
element to the state of minimum energy density Sq. These 
forces result from the elastic strain. The deformation 
energy is determined by the strain tensor, 



du, 
dxk 



duk 
dxi 



(15) 



and the elastic coefficients of the crust. For an ideal bcc lat- 
tice, there are two independen t elastic she ar moduli bn — 
\{cii — C12) and C44 (see, e.g.. lKitteilll986D . For a classica l 
Coulomb l attice, they were fi rst calculated bv lFuchsl (|l936l ) 
(see also lOgata et al.1 Il990t ): 611 = 0.0245nj^(Ze)Vrc, 
C44 = 0.1827n^(Ze)^/rc. Crustal matter is likely to be an 
isotropic bcc polycrystal, with a single effective shear mod- 
ulus /z. The shear strain contribution to the stress tensor is 
then 

1 



sh 



2^, 



E 



Ukk 



(16) 



A standard choice in the preceding studies of elas- 
tic effects in the neutron star structure and dynamics 
was u = CAA (iBavm fc Pines! Il97ll : IPandharipande et al.l 
119761 : iMc Dermott et al.lll988l and references therein). This 
choice significantly overestimates /i. A correct calculation, 
involvin g averaging oyer ro tations of Cartesian axes, was 
done bv lOgata et al.l (|l990t) . and gives 



li=\ {2bn + 3c44) - 0.1199 ''^^^''^' 

5 Tc 



(17) 



We rewrite this formula in a form suitable for calculations: 

^ = 7.816 X 102« {jf^y (I)' e^-g ' (18) 

where pia = p/10^^ g cm^'^, A^ is the total number of nu- 
cleons in the ion cell, Ac^s = Ac/10^, and Z is the num- 
ber of protons in an nucleus. This is the formula used in 



6 



J. L. Zdunik et al.: Deformation and crustal rigidity of rotating neutron stars 



the present paper. It should be stressed that this formula, 
strictly speaking, holds for the outer crust, where the ap- 
proximation of point-like nuclei is valid. However, for sim- 
plicity we also use Eq. in the inner crust, neglecting 
the effect of finite nuclear size and possible existence of 
non-spherical "pasta nuclei" (see the end fragment of Sect. 
SH). 



4.1. Ground-state crust 

This is a standard approximation assu med in the calcula- 
tions of the neutron-star structure (see lHaensel et aDl2007l 
and references therein). At each pressure, matter is assumed 
to realize the minimum of the enthalpy per nucleon. This 
is equivalent to the minimum of the energy per nucleon at 
every density; however, in contrast to the pressure that is 
continuous and strictly monotonous in the stellar interior, 
the density can suffer jumps. 

The ground-state approximation may be good for iso- 
lated neutron stars formed in core-collapse supernovae. 
Initial temperature is there ^ 10^° K and the matter is 
in nuclear equilibrium. In a few months, following neutron 
star birth, the temperature falls to 10^ K, so that layers of 
p ^ 10^ g cm~'^ are strongly degenerate. It is tacitly as- 
sumed that the matter cools, keeping nuclear equilibrium. 
As a result, the cooled crust reaches its ground state. Its 
structure and the EOS can be calculated assuming T — 0. 

In the inner crust, at densities approaching 10^'* g cm~'^, 
the nuclear radius becomes comparable to the distance be- 
tween nuclei. The effective nuclear charge becomes smaller 
than Z due to electrons present in nuclear volume, and 
moreover the nuclei can lose their spherica l shape, being 
repla ced by "nuclear pasta" (for a review, see lHaensel et al.l 
|2007| 1. However, the very existence of a "pasta nuclei" 
layer near the bott om of the crust is model-dependent 
(jHaensel et al■|[2"007^ . The pasta layers correspond to var- 
ious phases of liquid crystals and Eq. (flSj) is not valid 
for them. Elastic moduli of some pas ta phases were cal- 
culated by iPethick fc PotekhinI (|1998D (for a review, see 
iHaensel et al.ll2007f ). 

In the present paper, we limit ourselves to the crust 
models that do not exhibit pasta phases. W e mainly use the 
crust model of iDouchin fc Haensell (|200lh . in which spher- 
ical nuclei exist down to the bottom edge of the crust, so 
that Eq. (jlSp can be used as a reasonable approximation. 
In the upper panel of Fig. [71 we show plots of effective shear 
modulus /i versus density p for two rnodels of the ground- 
state crust (jPouchin fc Haensel l200lL iNegele fc Vautherinl 
Il973l the latter one for the inner crust only). Differences 
between two models are small. 



4.2. Accreted crust 

Neutron stars that have gone through a long stage of ac- 
cretion in the low-mass X-ray binaries have very different 
crusts from the ground-state one. Accretion could last for 
as long as ~ 10^ years, with typical accretion rates 10~^° — 
10-9 Mq y-i. After ~ 107(Af/10-9 M© y-^) years, the 
original crust, formed at neutron star birth, is replaced by 
an accreted one. During formation of an accreted crust, 
temperature in the stellar in terior with p ^ 10^ g cm"'^ is 
at most a few times 10^ K (jMiralda et al.|[T990f) . and the 
only reactions in compressed matter are those not associ- 
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Fig. 7. (Color online) Effective shear modulus p, Eq. p7|) . 
versus p, for different models of the cr ust. Upper panel: 



ground-state crust. Magenta line - Douchin fc Haensell 



green line - iNegele fc Vautherin (1973), only in- 
ner crust. Discontinuities in p have been smoothed via in- 
terpolation. Lower panel: blue line - an accreted crust 
mode l including pycnonuclear fusions (jHaensel fc ZduniS 
120081 ). The ashes of X-ray bursts consist of pure ^^Fe, 
the crust structure given in Table A. 3 of Appendix A of 
(jHaensel fc Zduni5l2008l ). Black line - ground-state crust, 
calculated using the same compressible liquid drop model 
of nuclei as the accreted crust. 



ated with overcoming the nuclear Coulomb barriers. Before 
the neutron drip point at p « 6 x 10^^ g cm^^, evolution 
is mediated by the electron captures. At higher density, 
neutron emission and absorption also play important role. 
Finally, at p ~ 10^^ — 10^^ g cm"^, the charge of proton clus- 
ters becomes so low that it has been speculated that their 
zero-point vibrations around the crustal la ttice s ites c ould 
possibly allow for the pycnonuclear fusion (jSatol . [19791 ). 

Accreted crust, produced in the scenario described 
above, has composition dramatically different from the 
ground-state one. In partic ular, nuclei i n an accreted crust 
have much sm aller A and Z (jSatol . [19791 : IHaensel fc Zdunikl . 
Il990al l2003t ). Therefore, an accreted crust is softer to an 
elastic deformation than is the ground-state one. To cor- 
rectly calculate the difference between shear moduli of the 
ground state and of an accreted crust, one has to use the 
same nuclear model. We a ssume the Macki e -Baym model 
of nuclei in de nse matter ([Mackie fc Bayrnl . 1977f) . in the 
version used in I Haensel fc Zdunikl (jl990aibl . l2003f) . 

Our results are presented in the lower panel of Fig. [71 
At a given p, nuclei in an accreted crust have lower values 
of Z and A than those in the ground-state crust. The val- 
ues of p{p) are correspondingly lower. One notices sharp 
drops and jumps of p of accreted crust. The drops are due 
to electron captures leading to Z — > Z — 2. In the in- 
ner crust, electron capture triggers emission of neutrons, 
so that both Z and A drop at the threshold density. At 
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some densities above 10^^ g cm^'^, Z may be low enough 
to allow for pycnonuclear fusion, resulting in a significant 
jump in Z and A, corresponding to a very large jump in 
/i. However, one has to stress uncertainties in the descrip- 
tion of zero point motion of proton clusters carrying with 
them bound neutrons, in the dense superfluid medium of 
unbound neutrons. Therefore, this part of the scenario of 
the formation of accreted crust is plagued by huge uncer- 
tainties, and even doubts, about whether the pycnonuclear 
fusion occurs at all. In contrast to the accreted crust, we get 
a smooth [i{p) for the ground-state crust, calculated using 
the same compressible liquid drop model of nuclei. Notice 
that our model has no shell correction term. 



5. Crust structure: an approximation 

In the simplest approximation, we consider a spherical neu- 
tron star built of an ideal fluid. The hydrostatic equilibrium 
is then determined by th e Tolman-Oppenheimer-Vo lkoff 
(TOV) equations (see, e.g.. IShapiro fc Teukolsk^ll983f) 



dP 
dr 



G{p + P/c^){m + A-KT^P/c^) 
r2 (1 - 2Gm/rc2) 



dm 
dr 



47rr p , 



(19) 



(20) 



where r is the circumferential radius and mir) the (gravi- 
tational) mass contained within the sphere of radius r. The 
boundary conditions are r7i(0) = and P{R) = 0, where R 
is the star radius. The star mass is M — m{R). The bottom 
of the crust is at r = i?b, so that the mass of the crust is 
Mcr = M — m(i?b). The pressure at the bottom of the crust 
is Pb =P(i?b). 

Within the crust, < P < Pb, R > r > and M - 
Mcr < m < M. We limit ourselves to neutron stars with 
Mcr < M, so that m^i M and 'iirr'^ P / mc^ < 1 within the 
crust. Moreover, within the crust P/c^ <^ p. Under these 



(21) 



conditions, Eq. (jl9p can be rewritten as 



dP 
dr 

or, equivalently, as 

dP _ 
P 



GMp 



r2 (1 - 2GM/rc^) 



-GM 



dr 



r2(l - 2GM/rc^) 



(22) 



Within the crust, i.e., for < P < Pb, we define a dimen- 
sionless function of pressure 



, , dP' 



(23) 



Notice that x{P) is determined solely by the EOS of the 
crust. Using Eq. ((^ . one obtains x versus r, 



1 -r„/P 
1 — To- /r 



(24) 



where rg = 2GM/c^ . Defining a = 1 — r^/R, we obtain r 
within the crust as a function of x. 



= rg/ (1 - ac-2x) . 



(25) 



We can now derive a useful formula for an integral of any 
function / over the crust volume. A proper volume of a 
spherical shell between r and r -I- dr is 

dV^inr^- ^^^^r^TTTT . (26) 



(l-rg/r)V2 



Using Eq. ((25|) we obtain 



r fdV^8nR'{l~afa^^^ f f- 



dr] 



arj 



2\A 



(27) 



where 77 = e~-^ and 77b = 77(Pb). Therefore, for a given 
function /(P) and for a given value of Pb, it is sufficient to 
calculate a function 



(28) 



and the value of the integral of / over the crust of any star 
having the radius and the mass of the star R and M can 
then be obtained using Eq. ((27)l . In our case f — p. 



0.5 
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- BNIHI 
GNHQmS 



0.4 



0.6 



0.8 



Fig. 8. Upper panel: the function Z^{a), Eq. ([28]), calcu- 
lated for the DH model of the ground-state crust, assuming 
p{P) — p{P), Eq. pT|) . Lower panel: relation between 
gravitational mass M and a — 1 — rg/ R, calculated for four 
EOSs of the neutron-star core. 



6. Calculation of B 

We assume a quadrupolar deformation of the star. Starting 
from a relaxed rotating configuration C, characterized by 
some Cl and e, we calculate the hydro-elastic equilibrium 
configuration C^'^ after a continuous, complete spin down 
to O = 0. The configuration C^'^ is characterized by > 0. 
Within our approximation, the hydro-elastic equilibrium at 
= implies 



6 = 



B = 



-A K bA. 



(29) 
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Classical analytic result, valid for a self-gravitatin g elas- 
tic sp here of constant p and /x, and volume V ^ is ()Lovd . 
119201) 

^7 

B^^^V,^ ,30) 

In this case the dis placement field nfr) is calculated analyti- 
cally (|Lovelll920D . lBavm fc Phlei (|l97lh used this displace- 
ment field to c alculate B = Bbp for a compressible solid 
crust (model o flBavm et al.lll97"lbl). ei iveloping a liquid core 
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with EOS of iBavm et alT (|l971a[) . iPandharipande et "all 
(|1976| ) performed calculations of B using t wo approxima- 
tions. First, they used the prescription of iBavm &: Pinea 
(|l97 lf). Second, they calculated B using a modification of 
Eq. ([50)1 obtained by replacing 57/50 by one, V by crust 
volume Vc, and by the average over the crust volume, JI. 
Therefore, their approximation reads as 



B w B 



int 



d^ r/^(r) 



(31) 



IPandharipande et alj (|1976[ ) found that -Bjnt is a rather 
good approxima tion of Bbp- 

As shown bv lCutler et al.l (|2003[ ). the use of u(r) of con- 
stant p and p solid sphere for a thin solid shell envelop- 
ing a fl uid core badly overestimates B. Using the analytic 
ujr) of Franco et al. I ((200 0') for an incompressible shell, 
ICutler et al. ( 2003[ ) derived, within the Cowling approxi- 
mation, a correct formula for a thin crust model of constant 
p and (U, 

12 



= 55^^^^ 



(32) 



The prefactor is smaller by 5.23 than in Eq. ([50]) . However, 
the use of Eq. (|32p with a simple substitution p — > Jl 
would still badly overestimat e the true va l ue of B, by a 
factor « 8 for a 1.4 Mq star (jCutler et all . l2003l ). This is 
due to additional cancellations of different contributions of 
elastic strain to e,. ( Cutler ct al., 2003). The calculations 
of ICutler et^ (I2003D were done using the model of the 
crust of iNegele fc Vautherinl (|1973D . combined with several 
E OSs of the core. Nu merical results, presented in Table 1 
of ICutler et all((200l . were obtained for a 1.4 Mq neutron 
star. 

We do not attempt to find B vi a solving the hydro - 
elastic equilibrium equations, like in ICutler et al.l (|2003D . 
Instead, we approximate -B by a generalization of Eq. (|3T|) 
in the form 

B = 0.027 / d^r p{r). (33) 

J crust 

where the prefactor accounts for the real distribution 

of strain stresses within the crust after the spin down. 

Assuming a reduction factor 0.027 « 1/42 as compared 

to the Baym-Pine s' estimate, we reproduce results of 
ICutler et all (|2003f ) quite well. 

7. AB(/9b) for the ground-state and accreted crust 

We denote by AB{pb) the contribution to B of the crust 
layer with density pb at the bottom. We use the ap- 
proximation given by Eqs. p3|) and (|27|) . Therefore, the 
value of AB(p) is determined by two functions, crust EOS 
and p{P), and by a single, star-structure parameter a = 
1 — Tg/i?. The dependence of AB{p) on the crust model 
is illustrated in Fig. jS] Three curves were obtained using 
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Fig. 9. (Color online) Function Ai?(pb) calculated for three 
models of neutron-star crust and assuming a = 0.65. This 
value of a is obtained for a 1.4 Mq s tar model based on 
the DH EOS for the crust and core (jPouchin fc Haensell 
i200ll ). Two MB models differ in the treatment of the pyc- 
nonuclear fusion. In both MB models. X-ray ashes are pure 
^^Fe. Green dotted line: artificial blocking of pycnonuclear 
fusion until 2^min = 4. Red dashed line: no artificial sup- 
pre ssion of pycnonuclear re actions. Table A. 3 of Appendix 
to (|Haensel fc Zdunikll2008f) . 



different crust models. The upper curve was bas ed on the 
ground-state model of iDouchin fc Haensell (j2001f ). The in- 
termediate curve was calculated ass t iming a model of ac- 
creted crust from lHaensel fc ZdunikI (|2008( ) (see caption to 
Fig. ini) . The values of p are lower than for the ground-state 
crust, making AB[p\,) lower, too. 

Pycnonuclear reaction rates at p > 10^^ g cm~'^ are 
the most uncertain element of an accreted crust model. 
An uncertainty factor of ^ 10^ is common for theoret - 
ical calculations of reaction rates (|Yakovlev et al.l [20061 ). 
To visualize the role of pycnonuclear fusion in determin- 
ing Ai?, we (artificially) suppressed these reactions for 
Z > Zmin — 4. Such modification has no significant effect 
on the total deep crustal heating in accreting neutron star 
(|Haensel fc Zdunikll2008[) . However, as seen in Fig.jSJ it has 
a strong effect on Ai?(pb), decreasing it by a factor of two 
at 2 X 10^'^ g cm~^ and by a factor of five at 10^"* g cm~^. 
By nearly doubling Z and decreasing the number of nuclei 
by two in each reaction shell, pycnonuclear reactions, have 
a very strong impact on the crustal rigidity. 



8. Scaling of B with M and R 

iJonesI (j2000|) derived, within Newtonian gravity, a simple 
scaling of /^.^^^j /LtdF with M and R. The starting point 
was an estimate of the crust thickness Ai?, in the limit 
of AR/R ^ 1. By balancing gravity and pressure at the 
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exact 

integral approximation 
thin crust approx. 




Fig. 10. (Color online) Parameter B versus stellar mass 
M calculated for the DH EOS (Table 1). Solid line: exact 
result obtained using Eq. (j33[) . Dashed line: approximation 
explained in Sect. [H Thi n do tted line: scaling with M and 
R, as derived bv lJoned (j2000t) . normalized here to the exact 
value of B at M = Mmax- 



bottom of the crust, one gets 

AR oc R^/M . 



(34) 



In the same approximation, the crust volume scales 

AV = A-kR^AR cx R^/M . (35) 

Defining a mean value of the shear modulus, ave raged over 
the crust volume, (M)crust' ^"^^ S^*^ ()Jonesll20"00l ) 



i?tc oc AV « RVM 



(36) 



where the index "tc" recalls that the formula is obtained in 
the limit of a thin crust. In what follows, we "normalize" 
the prefactor in the formula for Btc by requiring that the 
formula reproduces (for a fixed EOS for the core and the 
crust) the value obtained by an exact integration /^.^^gj /^dV^ 
(see the end of this section). 

In Sect.[5]we derive an approximate, but as we will show 
very precise, formula for /^.^.^^^ ^AV . In particular, we found 
that for a given n{P) and EOS of the crust, the relevant in- 
tegral over the crust scales in a well-defined manner with R 
and M, the latter scaling entering through the dependence 
on a = 1 — Tg/R. If we change M, R — > M', i?', we get, in 
our approximation 



B^B' = {R'lRfZ^{a')/Z^{a) 



(37) 



where 2fj_{a) is given by Eq. 

To be specific, let us consider the DH EOS of the crust 
and the core, for which we also calculated fi = fJ-{P)- We fix 
the prefactor in Btc by requiring Btc = B at M = Afmax, 
where the crust is the thinnest. In Fig.[TO]we compare B{M) 
obtained using three approximations. The method based on 



Eq. (P7)l is very precise, even for M « 0.3 M© when the 
crust is no longer thin, while for M > 0.5 Mq its precision 
is remarkably high. Approximation Btc is reasonably good 
only for M > Mq and provided the value of (m) crust 
properly chosen, which requires a "normalization" by an 
exact value of the integral at a selected value of M. For 
M < Mq, Btc diverges badly from the precise B{M), and 
the scaling given by Eq. ([5^ cannot be used. 



9. Numerical results for crustal rigidity b 

Our results are presented in Fig. [TT] In the upper panel we 
show our results for b{M) calculated using the DH model 
of the crust and core, for 0.25 Mq < M < Mmax(DH) = 
2.05 Mq. In the lower panel we have plotted b{M) for M > 
1 Mq , calculated using the DH model of the crust and four 
different EOSs of the core. 




-0.2 

log(M) [MJ 
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Fig. 11. Crustal rigidity b versus stellar mass. Upper 
panel: log-log plot, which shows results for the DH model of 
the ground-state crust and core for 0.2 M© < M < Afmax- 
Lower panel: zoomed segment of b{M) for M > 1 Mq 
(within a rectangle in the lower-right corner of the upper 
panel) for the DH model of the ground-state crust combined 
with several EOSs of the core. 



9.1. Nucleon cores 

As examples, we used the DH and APR EOSs for the 
core. The core EOS does not exhibit any strong soften- 
ing due to hyperonization or a phase transition, so b(M) 
is monotonous. It increases with decreasing M because of 
a monotonous decrease in A (Fig. |4]) and an increase in 
B due to an monotonous increase in the crust volume. 
Crustal rigidity is a sensitive function of M; the value of b 
at M = 1.0 Mq is ten times greater than at M = 2 Mq. 
For a "standard" mass 1.4 M© we get b ^ 2 x 10^^, which 
is within the range of 6 = 1.5 — 4.4 x lO"*", obtained by 
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ICutler et all (|2003[ ) for the FP S crust (fLorenz et al]ll99| ) 
and Negele-Vautherin crust (|Negele fc Vautherinl Il973f) . 
each com bined with two dif ferent EOSs of the core (see 
Table 1 of lCutler et aI1l2003h . 

The increase in b with decreasing M is so strong that 
6(0.3 Mq) w 506(1.4 Mq) and we were forced to use the log- 
arithmic scale to plot 6 in the whole considered mass range. 
Even at M = 0.3 Mq, the maximum elastic strain effect 
on neutron star structure, measured by 6 ~ 10~^, is a very 
small perturbation of the perfect-fluid model. However, one 
also has to notice that the applicability of the small defor- 
mation approximation to neutron stars with Al < 1 Mq, 
rotating at a few hundred Hz, can be questioned. 

Consider, for example, / = 641 Hz, corresponding to 
PSR B1937-f26. For illustration, wc use properties of the 
low-mass neutro n stars calculat ed for the DH EOS of the 
crust and core bv lHaensel eFaTI p00'2( ). At / = 641 Hz, the 
equatorial radius of a 0.6 Mq star is 17.8 km, to be com- 
pared with 12.6 km for a nonrotating star of the same mass. 
The polar flattening is characterized by rpoic/j'cq ~ 0.6 
(here r is the radial coordinate). The equatorial thickness 
of the crust is four times that at the poles. Actually, at 
/ = 641 Hz, stars with M < 0.6 Mq cannot exist in hy- 
drostatic equilibrium: they are destroyed due to the mass 
shedding from the equator. 

If DH EOS of the core is replaced by the APR one, the 
change in 6(M) is small. Namely, 6(M) is somewhat lower, 
due to a higher stiffness of the APR EOS. 

9.2. Hyperon cores 

A strong softening of the core EOS due to hyperoniza- 
tion has a dramatic effect on the high- A/ segment of b{M). 
Namely, the crustal rigidity is no longer monotonous, and 
after reaching a rather narrow maximum of nearly 10^^ at 
~ 1.5 Mq, it falls monotonously, down to about 10~^, at 
M = A'/max(BGNlHl) = 1.65 Mq. This behavior is easily 
understood in terms of a specific dependence of A on M 
(see Fig. [5]) , characterized by a deep and narrow minimum 
around ^ 1.5 Mq. A minimum in A{M) is reflected by a 
maximum in b{M). 

9.3. Exotic cores 

Exotic cores are represented by the GNHQm2 EOS, with a 
phase transition to quark matter via a mixed quark-hadron 
state. The softening effect implies a low value of Mmax, high 
compactness of "hybrid stars" and rapid shrinking of the 
crusts with M approaching Mmax- This results in a rapid 
drop in 6 near Mmax- The effect of a softening due to mixed 
state is seen in Fig. [5l It is represented by flattening of 
the A{M) curve GNHQm2, and it is more gradual than 
a violent softening by a hyperonization, as seen in curve 
BGNIHI in the same figure. 



10. Discussion and conclusion 

In the present paper we have studied the dependence of 
the A and B parameters of the Baym-Pines model of the 
hydro-elastic equilibrium of rotating neutron stars on the 
stellar mass and on the EOS. We considered seven different 
EOSs of the liquid core and several models of neutron star 
crust. 



Parameter A. It was calculated using high-precision 2-D 
simulations in General Relativity, so it can be considered 
as exact, in contrast to the previous estimates based on 
the incompressible fluid model. We developed a systematic 
procedure for calculating A{M). This procedure is based 
on linear expansions of global parameters of rotating neu- 
tron stars of a fixed baryon mass in the square of rotation 
frequency as measured by an observer at infinity. Three 
different coefficients of the quadratic term measure, respec- 
tively, the anisotropy of the inertia tensor (a/), oblate- 
ness of shape (as), and increase in the energy as mea- 
sured in the star's rest frame (ad). Combined, they give 
A = ad{a.s/ aiY ■ A precise calculation of A{M) was pos- 
sible for 0.2 Mq <, M <, 0.9 Mmax- Our exact values of 
A{M) are usually much higher than those given by the 
standard formula based on the incompressible fiuid model, 
used in the previous calculations. For M < 0.7 Mq the 
values of A{M) are determined by a low-density segment 
of the EOS (similar for all EOSs) and are given by a uni- 
versal formula of precision on the order of one percent. For 
M > 0.7 Mq, values of A{M) are EOS dependent and di- 
verge with increasing M . For EOSs without a high-density 
softening (only nucleons and leptons), the values A{M) in- 
crease monotonously with M. For EOSs with a high-density 
softening (due to hyperons or a phase transition), A{M) can 
have a strongly non- monotonous high-mass segment, where 
A can drop strongly with an increasing M. For example, 
softening by hyperons in the BGNIHI EOS makes A drop 
by a factor of about three for M increasing from 1.4 Mq 
to 1.5 Mq. Generally, the EOS dependence of A{M) for 
M « 1.5 Mq is very strong, and A can increase by a factor 
of five when going from hyperon-softened EOS to a stiff 
nucleonic EOS of the APR type. 

Crustal rigidity. We derived an approximate formula for 
the integral of the elastic strain over the crust volume, 
whose validity goes beyond a thin crust approximation. Our 
formula can also be used for low-mass neutron stars. We 
considered three different models of the crust, two for the 
ground-state crust and one for an accreted crust. The values 
of B for different models of the ground-state crust and the 
same core EOS show little model dependence. In contrast, 
the accreted crust is significantly less rigid (lower B) than 
the ground-state one. Results obtained for M — 1.4 Mq 
with the ground-state DH crust and nucleonic cores ar e 
within the range of values obtained bv lCutler et al.l (|2003l ): 
however, we found a strong dependence of 6 on the star's 
mass. In the range 1.0 Mq — Mmax, parameter 6 varies 
by an order of magnitude, and one notices a characteristic 
dependence of relation 6(M) on the EOS of the core. For 
nucleonic cores, b{M) is monotonic. However, hyperoniza- 
tion may lead to a strong maximum in b{M), associated 
with soft hyperon cores and a deep minimum in A{M). A 
phase transition via a mixed quark-hadron state implies a 
strong drop in 6. 

Pycnonuclear reactions and crustal rigidity. We have shown 
that pycnonuclear reactions, if they proceed, play a cru- 
cial role for rigidity of accreted crust, because every time 
they occur they nearly double Z. However, theoretical py- 
cnonuclear reaction rates are notorious for their huge un- 
certainties. There is also a basic difficulty in describing of 
aggregates of neutrons, bound to the proton clusters, bound 
in turn in the nuclear potential well created by neutrons, 
as "nuclei" immersed in a superfluid of dripped (unbound) 
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neutrons. The effect of dripped neutrons on the zero-point 
vibrations of such "nuclei" and on the fusion rate stiU needs 
to be calculated. We use a compressible liquid drop model 
of nuclei in the inner crust. Assuming that pycnonuclear fu- 
sions are suppressed for Z > Zj^in = 4, we get B, which is 
five times smaller than for a model without such an artificial 
suppression. This may be contrasted with deep crustal heat- 
ing in accreting neutron stars, where the total heat release 
per one accreted nucleon was found to b e insensitive to the 
supp ression of pycnonuclear reactions (|Haensel fc ZduniS 
I2008D . 

This paper is first in a series devoted to the inter- 
play of rotation and shear strain in neutron stars. The 
forthcoming papers will deal with hydro-elastic equilibrium 
of rotating neutron stars, including the case of hypotheti- 
cal solid neutron-star cores, predicted by some theories of 
dense matter. Elastic strain can support non-axial mass 
distributions of rotating neutron stars, therefore becoming 
sources of continuous, strictly periodic gravitational waves. 
Such models of gravitational wave sources are of particular 
interest in the expected era of gravitational-wave astron- 
omy, which could be (hopefully) opened by LIGO, VIRGO, 
and other gravitational-wave observatories in the coming 
decade. 
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